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ABSTRACT 

Several X-ray observations have recently revealed the presence of reverberation time delays between 
spectral components in AGN. Most of the observed lags are between the power-law Comptonization 
component, seen directly, and the soft excess produced by reflection in the vicinity of the black hole. 
NGC 4151 was the first object to show these lags in the iron K band. Here, we report the discovery 
of reverberation lags in the Fe K band in two other sources: MCG-5-23-16 and NGC 7314. In both 
objects, the 6 — 7 keV band, where the Fe Ka line peaks, lags the bands at lower and higher energies 
with a time delay of ~ 1 kilo-seconds. These lags are unlikely to be due to the narrow Fe Ka line. They 
are fully consistent with reverberation of the relativistically-broadened iron Ka line. The measured 
lags, their time-scale and spectral modeling, indicate that most of the radiation is emitted at ~5 and 
24 gravitational radii for MCG-5-23-16 and NGC 7314 respectively. 
Subject headings: AGN etc 



1. INTRODUCTION 

Observations of active galactic nuclei (AGN) have 
shown that most of the X-ray radiation is emitted 
very close to the central supe r-massive black hole . 
Short time scale variability (e.g. Vaughan et al. 2011), 
relativistic spectral distor t ions (Tanaka et al. 1995; 
i Revnqlds and Nowakl 120031 : iMillerl 120071: IFabian et all 
2012b^ X-ray eclipses ([Risaliti et al.l 12007ft . as 
well as gravitatio nal micro-lensing measurements 
(jChartas et al.l l2012f h all indicate that emission origi- 
nates in a compact region a few gravitational radii from 
the event horizon. 

The primary X-ray radiation is thought to be pro- 
duced in a hot electron coro na by Compton up-scatterin g 
lower energy disk photons (jHaardt and Marasc hi 1991). 
This radiation illuminates the su r round i ng, relatively 
cold, matter (jGuilbert and Reesl 119881 : IFabian et al.l 
1989 ) giving rise to c haracteristic reflection spectra 
(R oss and Fabianl 119931 ). The variability of the pri- 
mary X-ray source in this environment naturally leads 
to time delay s comparable to the light-crossing time to 
the reflector (Re ynolds et al.| ["l999). The first direct evi- 
dence for such delays was measured in the Narrow Line 
Seyfert 1 (NLS1) galaxy 1H0707-495 (|Fabian et al.ll200l 
Zogh bi et all [201(1 . It was observed as a ~ 30 sec- 
onds delay between the primary X-ray source and the 
strong iron L line produced by reflection. Several sub- 
sequent studies have shown that such lags are present 
in many other Seyfert ga l axies p ossessing a soft excess 
(lEmmanoulopoulos et all 120111: [d e Marc o et all 120111: 
IZogh bi and Fabian 2011; Kara et a l. 2013b: [Fabian et"aU 
l2012aHCackett et al.ll2012l ). All these delays are compa- 
rable to the light-crossing time at a few gravitational 
radii (r g ) from the black hole, and there is even evidence 
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that the amplitude and f requency of the lag co rrelates 
with the black hole mass (De Ma rco et al.ll2012ft . 

A key property in all those objects is their strong soft 
excess, a characteristic feature in the NLS1 class. The 
measured time delays between the direct power-law com- 
ponent, that dominates the 2 — 4 keV band, and the soft 
excess (< 1 keV) are naturally explained if the latter is 
due to a reflection process. In this case, a combination of 
emission lines and bremsstrahlung from the surface lay- 
ers of an ionized disk are smoothed and distorted by rel- 
ativistic effects. Interpreting the measured lags as light- 
crossing effects in t he immediate vicinit y of the black hole 
was questioned bv lMiller et al.l ((20 10a) who argued that 
these delays could be an artifact of a much larger sys- 
tem (hundreds of r g ). Although such an interpretation 
fails when a full set of obse rved timing and spe ctroscopic 
properties are considered ([Zoghbi et al.l 1201 If ), it partly 
stands on the ambiguity of interpreting the soft excess it- 
self. The spectrum of the observed soft excess tend to be 
very smooth masking out any direct detection of spectral 
emission features. Detecting reverberation of the broad 
iron K- alpha line at ~ 6 keV would remove such ambigu- 
ity and clearly demonstrate an accretion disk origin for 
the time lags. 

Lags between energy bands containing the iron-K line 
and softer energy bands have been k nown about for 
some time, both in Gala ctic black holes (Mivamoto et al.l 
1988 :lKotov et al.ll2001ft and AGN (jMcHardv et alJl200ft 
lArevalo et al.l [20081 ). They are seen at low frequencies 
(typically around 10 -4 Hz for a 10 6 M black hole, 
and lower freque ncies for higher ma s ses or larger emis- 
sion radius, e.g . IZoghbi et al.l 120111 : iKara et all l2013at 
iDe Marco et al.l 1201 2ft and show an almost featureless 
energy spectrum (with At ~ logE); which, in analogy 
with Galactic black holes (where the lags extend well 
above 10 keV), leads most authors to associate them 
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with continuum process it self (e.g. Comptonization or 
propagating fluctuati ons. iKotov et al. 20 Oil [Poutanenl 
120011 : iMcHardv et al.l 12007 : Zogh bi et al.ll2Qlih . On the 
other hand. iMiller et al.l (|2Q10al lbl) attribute such low- 
frequency lags to reverberation associated with the con- 
tinuum scattering in a hypothesized extended structure. 
A major recent development was the discovery that the 
lags in NGC 4151 which has an energy dependence 
that traces th e shap e of a relativistically-broad iron line 
(Zog hbi et all 120121 ). This is a clear and dramatic con- 
firmation of a central prediction of the relativistic disk- 
reflection paradigm. The observed lag in this case is 
an extension of the s oft reverberation lag s seen in many 
objects (jFabian et al.l 120091 : IZoghbi et aLll2010l and later 
work), and it is separate from the low- frequency feature- 
less continuum lag. 

In addition, the 'line' seen in the lag-energy spectrum 
of NGC 4151 has a stronger blue horn at long time-scales 
and a prominent red wing at short time-scales, in a strik- 
ing match to the expectation of a lin e emitted from a dis k 
and distorted by relativistic effects (Fabia net al.l fl989). 
Detection of Fe K l ags has also been seen in 1H0707-495 
(iKara et all l2Q13bl ) and IRAS 13224-3809 ([Kara et al.l 
l2013bl ). Here, we report on the detection of similar lags 
in two Seyfert galaxies MCG-5-23-16 and NGC 7314 us- 
ing XMM- Newton observations. The two objects were 
selected because of their brightness and high variability 
in the 2 — 10 keV band. Their selection came during the 
search for the best reverber ation targets for the proposed 
X-ray mission GRAVITAS dNandra et al.ll2Q12h . 

MCG-5-23-16 is a Seyfert 1.9 Galaxy (z = 0.0085) with 
a typical 2-10 keV flux of ~ 8 x 10~ n ergs cm" 2 s -1 
and a mass of 5 x 10 7 M (jWandel and Mushotzkvl 
1986). Its spectrum resembles a classical Compton- 
thin Seyfert 2 galaxy with a column that does not af- 
fect the spectrum above 3 keV significantly. The spec- 
trum below 1 keV contains a combination of emission 
from scattered continuum photons and distant photoion- 
ized gas. Above 2 keV, the spectrum shows both a 
narrow (EW-60 eV) and broad (EW~ 50 - 200 eV) 
iron Ka l ines along with a s t rong Compton hump above 
10 keV dWeaver et al.l 119971: iMattson and Weaver! 120041: 
iBraito et al.l 120061 : iReeves et al.ll2007D . 

NGC 7314 is also classified as a Seyfert 1.9 galaxy 
z = 0.0048) with an estimated mass of 5 x 10 6 M 
Schulz et al.l Il994f ) . It shows strong variability on all 
observed time scales. It i s thought to be a type 2 counter - 
part to the NLS1 class (Dewangan and Griffiths] l2QQ5h . 
The source is possibly seen through a weak warm ab- 
sorber. The iron K band also shows both a narrow and 
broad com ponents in observations with ASCA, Chandr a 
and XMM (jYaqoob et al.lll996Ll2003l : lEbrero . et aLl2011h . 
It has also been shown that the narrow and broad com- 
ponents of t he line respond differently to variations in the 
continuum (jYaqoob et al.ll2QQ3h . Both objects show hard 
and absorbed spectra. We concentrate in our analysis on 
energies > 2 keV to avoid any spectral complexities due 
to either emission from scattered and distant continuum 
photons and photoionized gas, or from warm absorption. 

The paper is organized as follows: Section [2] de- 
scribes the data used in the analysis. The results for 
MCG-5-23-16 and NGC 7314 are presented in Sections 
[3] and |4] respectively, which include both lag and spectral 
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Figure 1. The light curve of MCG-5-23-16 in the 2-10 energy 
band. It shows the high count rate and strong variability. 

modeling. The interpretation and implications of the re- 
sults are presented in Section [5] 

2. DATA ANALYSIS 

The results presented here are obtained using archival 
XMM observations of MCG-5-23-16 and NGC 7314. 
MCG-5-23-16 was observed three times (Obs ids: 
0112830301, 0112830401 and 0302850201). The first had 
a strong particle background and is not used. The sec- 
ond observation is short (20 ks after high background 
filtering) and not very useful for timing analysis. Only 
the third long observation was used for timing (expo- 
sure: 130 ks) after removing the last ~ 20 ks affected by 
strong particle background in the detector. This obser- 
vation was taken in the large window mode and had some 
pileup, with a deviation from the model by 1.6 and 3.7 
percent for the single and double patterns respectively 
(as measured by epatplot in SAS) . Pile- up is not an 
issue for timing analysis as show by Zogh bi and Fabianl 
(2011), who performed detailed simulations on the ef- 
fect of pile- up on lag measurements. The central regions 
were exercised in the spectral modeling to remove the 
most affected regions in the psf. NGC 7314 was ob- 
served only once (obs id: 0111790101; exposure: 44 ks) 
in a small window mode with an additional one off-axis 
observation (obs id: 0311190101; exposure 83 ks) and 
both observations were used. The ODF files were reduced 
with the standard pipeline using SAS 12.0.1. Source and 
background light curves were extracted from circular re- 
gions of 50 arcsec radii from the pn detector. The light 
curves were then corrected for instrumental effects using 
epiclccorr. Spectra and responses were extracted fol- 
lowing the standard procedure in the XMM ABC guide. 

3. MCG-5-23-16 

3.1. Time lags 

Fig. [1] shows the 2 — 10 keV light curve. The object is 
clearly very bright and variable. To study time lags, light 
curves were constructed in eight energy bins between 2 
and 10 keV with 1 keV bin width. We used standard 
Fourier tec hniques to ca l culate the frequency-dependent 
time lags (Nowa ket al.l [1999) and found that the lag 
depends on both frequency and on the energy bands 
used. As an example Fig. [2] (Top) shows the frequency- 
dependent lag between the 4 — 5 and 6 — 7 keV energy 
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Figure 2. The lag (top) and coherence (bottom) as a function of 
Fourier frequency between 4 — 5 and 6 — 7 keV energy bands for 
MCG-5-23-16. Positive lags in the top panel indicate hard lags, 
where the harder band (6 — 7 keV, where the reflected iron K line 
peaks) lags the softer band (4 — 5 keV, dominated by the direct 
emission). There is a significant hard lag at frequencies < 10 -4 
Hz. The measured coherence remains high at those frequencies. 

bands. We use the standard sign convention, where pos- 
itive values indicate a hard lag (i.e. the harder band is 
delayed with respect to the softer band). It is clear that 
the 6 — 7 keV band lags behind the 4 — 5 keV band below 
10 -4 Hz. The bottom pane l of Fig. [2] shows the mea - 
sured coherence function 7 2 ( Vaug han and Nowak 1997). 
This is a measure of the fraction of one light curve that 
can be predicted from the other. It can be seen that the 
coherence is high and the light curves match each others 
variability. 

In order to study the energy-dependence of the lag, we 
select the frequency band where there is a positive lag 
in Fig. [2]-top, i.e. / < 10 -4 Hz, and calculate the lag 
for the different energy bands with respect to a reference 
band. The reference can be any band, or can be the total 
2 — 10 keV band itself. Here, for each energy band, the 
reference band is taken to be the whole 2 — 10 keV band 
excluding the current band. Using the whole band as a 
reference maximizes the signal to noise, while removing 
the band of interest from the r eference ensures th at the 
noise remains uncorrelated (see lZoghbi et al.ll2011f ). The 
result is a plot of the lag at frequencies < 10 -4 Hz as a 
function of energy, and it is shown in Fig. [3l 

Fig. [3] clearly shows that the lag traces the shape of 
a line that peaks in the 6 — 7 keV band, where the iron 
Ka line peaks. This is similar to the shape reported 
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Figure 3. The lag as a function of energy for frequencies < 10 -4 
for MCG-5-23-16. Energy bins have a width of 1 keV. The lags are 
measured with respect to a reference band taken to be the whole 
2 — 10 keV band (see text for details). 



for NGC 4151 (jZoghbi et al.ll2Q12[ ). but appears to be 
less broad. Note that the zero point in Fig. [3] (and the 
rest of the lag-energy plots in this work) depends on the 
chosen reference band, but for clarity, we have shifted the 
plot vertically so that the zero point is the point with the 
most negative lag (the last point in this case). Note also 
that the average value of the positive lag at frequencies 
< 10 -4 Hz in Fig. [2]-top is ~ 0.7 ks, which is the value 
of the vertical difference between the 4 — 5 and 6 — 7 
keV points in Fig. [3l To look for frequency-dependent 
change in the feature at ~ 6 keV, we also produced the 
same plot for the lowest frequencies (< 2 x 10 -5 ). Fig. 
[2]-top indicate that the lags becomes larger as we move 
to lower frequencies, and this is indeed what we see in 
Fig. H 
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Figure 4. The lag as a function of energy for frequencies < 10 — 4 
(red; continuous line) and < 2 x 10~^ (blue; dashed line) for 
MCG-5-23-16. The probability that the two lines have different 
peaks is 85%. 
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Table 1 

Best fit parameters. 



Parameter 


MCG-5-23-16 


NGC 7314 


r 

Fe abundance (Fixed) 


i kq -L n nfi 

-L.Oc/ ZC U.UU 

i 70 +o.03 

± -' U -0.13 
D.D_ 4 3 
1 


n aa A- n 1 ^ 
1.66 ±0.05 
243| 
1 


emm. index q 


2.2 ±0.5 


o 1+0.5 
Z - 1 -1.4 


Inclination angle(°) 


38i 9 8 


9 +io 


Gaussian Line 
Fe b road Eq. width (eV) 
X 2 /cLo.f 


93 ±13 
488/381 


6.92tn 9 6 keV 
82±g 
513/473 



The figure shows the lag for both frequency bands 
< 10 -4 Hz (red continuous line) and < 2 x 10 -5 Hz 
(blue dashed line). The lag magnitude increases at lower 
frequencies. There are even hints of a possible change 
of shape of the broad line, although not with a high sig- 
nificance. The probability of having different peak ener- 
gies is rsj 85% when fitted with gaussians with different 
widths. 

3.2. Spectral modeling 

To interpret the shape of the lag-energy plot, we turn 
to spectroscopy. The spectrum of MCG-5-23-16 above 2 
keV has been studied extensively. It was established that 
it has a broad iron lin e component, but n ot as extreme 
as some other objects. iReeves et al I ((2007) for example, 
using Suzaku and XMM data, measure a line width of 62 
eV, and an inner accretion disk radius of 37j^o r g when 
fitted with a diskline model (Fabia n~et al.| [l989). If a full 
reflection model is used, the inner radius is 26^g 5 r g if 
the emissivity index is constrained to be q = 3 (emissivity 
oc r _gf ), and goes to ~ 6r g if q is allowed to vary in the 
fit. Our fit for the XMM data gives 90% confidence upper 
limit on the inner radius of 18 r g and a best fit of 5.6 r g 
and a relatively flat emissivity index q = 1.8 ± 0.2. We 
used XMM spectrum only to remain consistent with the 
timing analysis. 

Fig. [5] shows the model used to fit the data between 
2 — 10 keV, and the best fitting parameters are presented 
in table [U There are two reflection co mponents modeled 
using the reflection table reflionx (Ro ss and Fabianl 
2005): The outer component has a characteristic narrow 
iron line and originates at large distances from the black 
hole (molecular torus or broad line region) and is not 
expected to produce any variability at t he tim e-scale of 
~ 10 ks probed by the lag analysis in sec. 13.11 The inner 
reflection component, modeled with the reflection table 
reflionx broadened by a relativistic kernel kdblur, is 
more likely to be responsible for the lags by responding 
to variability in the direct power-law component. 

Because both power-law and reflection contribute to 
every band between 2 — 10 keV, and if we assume, to 
a first order, the delay between the two components to 
be energy- independent, then each point in the lag plot in 
Fig. [3] gives a measure of the reflection fraction at that 
energy (variable reflection/power-law) multiplied by the 
intrinsic lag between the power-law and the reflection 
components (see further discussion in sec. 0) 




Energy (keV) 

Figure 5. The best fitting spectral model for MCG-5-23-16. 
The model includes a power-law component (blue dash-dot line), 
outer reflection (green dashed line) modeled with reflionx and a 
relativistically-broadened inner reflection (red dotted line), mod- 
eled with kdblur*ref lionx 

4. NGC 7314 

4.1. Time lags 

Fig. [6]-top shows the frequency-dependent lags for 
NGC 7314 between the the 4 — 5 and 6 — 7 bands. Similar 
to MCG-5-23-16, there is a hard positive lag, where the 
peak of the iron K line lags the 4 — 5 energy band, with 
a lag that depends on frequency. Fig. [6]- bott om shows 
the measured coherence function (see sec. 13. ip . The two 
bands have a relatively high coherence of ~ 0.8 up to 
- 5 x 10~ 4 Hz. 

Similar to the analysis in Sec. 13. 1\ we calculate the 
energy-dependent lags and show them in the Fig. [71 
there is a clear jump in the lag at ~ 5.5 keV. Its sig- 
nificance is more than 99%. The line profile traced by 
the la g is narrower than that of MCG-5-23-16 and NGC 
4151 (jZoghbi et al.ll2012[ ). We also tested for frequency- 
dependence of the profile. The narrow feature at 6 — 7 
keV does not change shape if we select frequencies say 
< 5 x 10 -5 Hz. There is only a change in the magnitude 
of the lag, which can already be seen in Fig. [6]-top. The 
constant shape with frequency is probably due to the 
sharpness of the jump ( or equivalently, the narrowness 
of the line ) making it difficult to detect a shape change. 

4.2. Spectral modeling 

We also performed spectral fitting of the XMM obser- 
vations. Previous spectral decomposition based on vari- 
ability showed that the sp ectrum at ^ 6 keV contains 
a narrow peak and wings (|Yaqoob et a l. 1996), indicat- 
ing that it possibly contains a broad comp onent. Using 
the hi gh resolution capabilities of Chandra. lYaqoob et al l 
(|2003f ) showed that the spectrum also contains emission 
from ions up to Fe xvil at ~ 6.4 keV, as well as transient 
emission features from high ionization species Fe xxv 
and Fe XX VI Ly<x All these lines appear to be variable 
and respond to continuum changes on in less than 12 ks. 
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Figure 6. The lag (top) and coherence (bottom) as a function of 
Fourier frequency between 4 — 5 and 6 — 7 keV energy bands for 
NGC 7314. 
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Figure 7. Energy-dependent lag for NGC 7314 at frequency of 
(1 - 30) x 1CT 5 Hz. The whole 2-10 keV band is used as a 
reference. 



Although it would be interesting to attribute the lags in 
Fig. [7] to a particular line or ionization state, the en- 
ergy resolution of the lags permitted by the current data 
(without significantly affecting the signal) cannot achieve 
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Figure 8. The best fitting spectral model for NGC 7314. The 
model includes a power-law component (blue dash-dot line), outer 
reflection (green dashed line) modeled with ref lionx (setting the 
ionization parameter to the lowest possible value of 1) and an inner 
reflection (red dotted line), modeled with laor and a narrow Fe 
xxvi line. 

that, and remains to be explored with future data. 

The spectrum from the XMM observation was ana- 
lyzed. There is a clear line at ~ 6.4 keV. If fitted with a 
gaussian, it has an energy of 6.43 ±0.02 keV and a width 
of 126 ± 35 eV, slightly resolved. To model it properly, 
we included the partially ionized reflection table model 
ref lionx in addition to an absorbed power-law. The 
residuals show the presence of a line at ~ 6.9 keV, which 
was modeled using a narrow gaussian . This could be 
due to Fe xxvi (see lEbrero et al.ll2011l for full spectrum 
analysis). There were also residuals at the wings of the 
6.4 keV line that are not accounted for by the cold re- 
flection model. This could be due to a broad component 
([Yaqoob et al.|[l996[ ) and we model the m wit h a gaussian 
or a broad relativisitc line laor (jLaorlll99lh . 

Although the broad component is not very broad, its 
presence is very significant (> 99.999% confidence when 
fitted with a gaussian). The inner radius of the emit- 
ting region from the laor model is not very constrained, 
but it extends down to 13 r g with an optimum value 
of 24 r gi with nearly face-on viewing angle of = 9° 
and a fla t emissivity o f q = 2, similar to the results 
found by lYaqoob et al.l (|2003f ) using Chandra observa- 
tions and different arguments. The best fitting model is 
shown in Fig. [8j and the best fitting parameters are pre- 
sented in table [U It should be noted that, although these 
are the optimum parameters, their uncertainties are rel- 
atively large. If more sophisticated models are used (e.g. 
kdblur*ref lionx), we obtain parameters that are simi- 
lar to those obtained from using the laor model. 

Although it is not our aim to do detailed spectral mod- 
eling, it is very clear that the spectroscopy is very con- 
sistent with lag measurements made in Sec. 14.11 (Fig. 
[7j), where there is a narrow line-like feature in the lag- 
energy plot. As we discussed in Sec. 13. 1| the lag-energy 
plot, under simple assumptions, gives a measure of the 
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reflection fraction as a function of energy (the variable 
reflection component responding to the direct power-law 
variations), and a reflection fraction that has a line peak- 
ing at ~ 6.4 keV is very consistent with the spectral 
modeling. 

5. DISCUSSION 

We have reported on the discovery of X-ray reverbera- 
tion in two bright AGN in the iron K band. The analysis 
in Sections 13.11 and 14.11 present model- independent mea- 
surements of time lags between the 6 — 7 band where the 
Fe Ka line peaks, and other bands, in a sense that this 
band always lags bands above and below. It is very ap- 
parent that these lags are linked to the emission in the Fe 
Ka band and very likely are due to a reflection process. 

The measurement of lag as a function of energy pro- 
vides a way of measuring the reflection fraction as a func- 
tion of energy in a completely independent way from the 
standard spectrum. This provides an invaluable tool to 
break the degeneracy often encountered in spectroscopic 
modeling, where it is often difficult to identify an un- 
derlying continuum due to absorption and emission com- 
plexities. To see how the lag-energy plot is meas uring the 
reflec tion fraction, consider the following (see iPoutanenl 
l200ll : iKotov et al.ll2Q0lh : Let S d (E,t) and S Te i(E,t) be 
the direct and reflected time-dependent spectra respec- 
tively, so that the total observed spectrum is: 



S(E,t) =S d (E,t) + S ie{ (E,t) 



(1) 



Assuming the direct component changes only in intensity 
and not in shape, we can write: 

S d (E,t) = D(E)A(t) (2) 

and the reflection can then be written as: 

S iel (E, t) = fR(E) J A(t - t)T(t)cIt (3) 

where D{E), R(E) are the direct and reflected spectra re- 
spectively, A(t) is the normalization of the direct compo- 
nent, T(r) is the transfer function of the reflected compo- 
nent and / is a constant that contain information about 
the geometry of the system. We have assumed again 
that R(E) only changes in amplitude as it varies. The 
frequency-resolved phase lag as a function of energy is 
then: 



tan(cP(E,f)) 



Im[S*(E,f)D(E,f)} 
Re[S*(E,f)D(E,f)} 



(4) 



where A denotes the Fourier transform and * denotes the 
complex conjugate. Re and Im indicate real and imagi- 
nary parts of a complex quantity respectively. This, for 
small reflection fractions then reduces to: 

tan(4>(EJ)) = f^Im[f(f)] (5) 

Therefore for a single frequency band, the lag as a func- 
tion of energy ( e.g. Fig. H]and[7j) gives a measure of the 
reflection fraction as a function of energy. 

Could these lags be due to the narrow Fe Ka compo- 
nent? Narrow Fe Ka lines are very common in AGN 
spectra. Although their origin is still not very clear 
(molecular torus, the broad line region or the outer parts 



of the disk itself, e.g. Nandra 2006), it is known that 
the narrow component of the line is not very variable 
generally, and certainly not on the time-scales probed 
by a typical single observ ation (~ 100 ks) such as 
those presented here (e.g. IVaughan and Fabianl 120041 : 
iBhavani and Nandral l2Q10h . For the data used in this 
work, we calculated the RMS spectra to assess the level 
of variability as a function of energy. We found that 
there is a drop in the RMS at ~ 6 keV. The width of 
the trough cannot however be constrained by the data, 
and cannot therefore be used to confirm or rule out the 
origin of variability in the narrow component. 

There are however other lines of evidence against the 
observed lags being due to the narrow Fe Ka line. The 
lags p resented here, and that of NGC 4151 ( Zo ghbi et al.l 
2012) appear in energy-dependence to be a continua- 
tion of the lags seen in the soft excess in many Nar- 
row Line Seyfert 1 galaxies. This is particularly ap- 
parent in the lag-energy dependenc e of 1H0707-495 
(iZoghbi e t al. 2011; Kara etaLl l2Q13bh , RE J1034+396 
(jZoghbi and F abian 2011) and IRAS 13224 (Kara et al.). 
Also, as in the case of NGC 4151, the observed peak in 
the lag-energy spectrum is less than 6.4 keV. If it was 
produced by the narrow component, the peak would be 
at the rest frame energy of the Fe line. Furthermore, 
the fact that the lag-energy shape depends on time-scale 
(i.e. temporal frequency) argues again against an origin 
in the narrow component. 

Now we consider the magnitude of the lags. The light- 
crossing time around a black hole at a radius r (in units 
of r g = GM/c 2 ) is bM^r seconds, where Mq is the 
black hole mass in units of 10 6 M . For M 6 = [50,5], 
where the first number is for MCG-5-23-16 and second 
for NGC 7314, the measured inner radii of emission in- 
ferred from the spectral fitting of r = [5.6,24] (in r g ) 
give a light crossing time of [1.4,0.6] ks, within a factor 
of 2 of the measured delays in Fig|3] and [3 It should be 
noted however that associating the measured lag directly 
with a particular radius is not straight forward. There 
are many factors that need to be considered. First, the 
measured lags in Fig. [4] and [71 include a dilution factor 
(equation [5]) that includes the reflection fraction which 
makes the measured lags between bands smaller than 
the intrinsic lag between the direct and reflected compo- 
nents. On the other hand, the effects of geometry and 
gravitational Shapiro delays make the radii inferred from 
light-crossing estimates larger than they are. For exam- 
ple, if the illuminating source is high above the disk, 
which could be the case for NGC 7314, then the delay is 
larger than the light-crossing time at the radius of emis- 
sion inferred from spectral fitting. Accounting for these 
effects requires full modeling that is beyond t he scope of 
this work (see e.g. IWilkins and Fabianl 12012} ) and would 
require more data to study the lag-energy in finer energy 
bins and frequency bands. 

Furthermore, as Fig. [2] and [6] show, the lag is 
frequency-dependent and a single measured lag at a sin- 
gle frequency cannot be trivially linked to a particular 
radius. It is however tempting to broadly associate small 
lags at small time-scales with small radii and longer lags 
with larger radii, The reason is that the competing effects 
of dilution vs geometry and gravitat ional delays roughly 
cancel out (though not exactly, see IWilkins and Fabianl 
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2012). Therefore, the increase in lag magnitude seen in 
Fig. |4]with decreasing Fourier frequency is likely caused 
by small time-scale variations at small radii (and hence 
small lags) being filtered out. The effect is similar to, and 
strongly supported by, the time-scale-dep endent shape 
of th e lag-energy spectrum of NGC 4151 (Zo ghbi et al.l 
12012ft , This highlights the power of these plots that mea- 
sure the reflection fraction independently of the classical 
spectrum, allowing us to use our knowledge of how the 
reflection fraction changes with time-scales to interpret 
the lags and the spectra. 

It is also interesting to examine the time-scale axis of 
the variability (i.e. Fourier frequency axis). Although 
the variability in black holes is well established by ob- 
servations, it is not clear what physical time-scale in the 
accretion disk it represents. For example, if we consider 
the viscous time-scales t v = 5r 1 - 5 a _1 (i^/ 'R)~ 2 Mq sec- 
onds, where r and Mq are again in units of r g and 10 6 
respectively, a nd a is the disk viscosity parameter 
(jFrank et al.l [2002). then the measured radii from spec- 
tral fitting of [5.6,24] r g give a viscous frequency l/t v 
of [7 x 1(T 8 ,8 x 10~ 8 ] Hz, using a = 0.1 and (H/R) = 
0.05. The measured lags clearly span higher frequen- 
cies (shorter time-scales) than this. If instead thermal 
time-scales are considered, where t t h = (H/R) 2 t Vl then 
the measured radii give frequencies of ~ 3 x 10 _5 Hz in 
both cases, which falls exactly in the observed range over 
which the reverberation lag is measured. Thermal time- 
scales could dominate if variability in density and/or ion- 
ization are important. Thermal time-scale seems to also 
be most relevant in the case of accreting white dwarfs 
too (e.g. iScaringi et al.ll2012l for some recent results). 

In summary, we have presented lag measurements in 
the Fe K band that show the peak of the line 6 — 7 keV 
band lagging bands either side of it. These simple mea- 
surements are model-independent, and to interpret them, 
we showed that they are consistent with a broad iron line 
responding to variations in an illuminating continuum. 
Most of the emission originates at 5 r g for MCG-5-23-16 
and 24 r g for NGC 7314. These two objects are selected 
for their brightness and high variability, two other vari- 
able N LS1 appear to show similar iron K lags (jKara et al.l 
l2013bl lal). which is a possible indication that these lags 
are common among variable objects. The question of 
how common are they among radio quiet AGN in gen- 
eral remains however to be explored. 
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